#include<iostream> 
#include<cylrrs.hpp>
//    Компиляция
//   g++  -Wall -I/usr/include -I/usr/local/include -I. -c cylrrs.cpp
//   g++  -I. -c main.cpp
//   g++  -L/usr/local/lib -o work cylrrs.o main.o -lgsl -lgslcblas -lm

int main(){
   cylrrs::init(3.76, .16, -60, 20000.0, 1.0/3, .25);
   //cylrrs::paraminfo();

   // Пример вывода компоненты $\sigma_\theta$ по слою
   // Вывод в единицах кГ/мм^2 Для перевода в МПа  необходимо умножить на 9.8
   double a = cylrrs::getradius();
   double dh = 0.01;
   for (int i=0; i<=20; i++){
      std::cout << i*dh << "  " << cylrrs::getsigmatheta(a-i*dh) << "\n";
   };
}
